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ABSTRACT 

In this paper, we examine the phase diagram of quenched QCD with 
two flavors of Wilson fermions, proposing the following microscopic pic- 
ture. The super-critical regions inside and outside the Aoki phase are 
characterized by the existence of a density of near- zero modes of the (her- 
mitian) Wilson-Dirac operator, and thus by a non-vanishing pion con- 
densate. Inside the Aoki phase, this density is built up from extended 
near-zero modes, while outside the Aoki phase, there is a non-vanishing 
density of exponentially localized near-zero modes, which occur in "excep- 
tional" gauge-field configurations. Nevertheless, no Goldstone excitations 
appear outside the Aoki phase, and the existence of Goldstone excitations 
may therefore be used to define the Aoki phase in both the quenched 
and unquenched theories. We show that the density of localized near-zero 
modes gives rise to a divergent pion two-point function, thus providing 
an alternative mechanism for satisfying the relevant Ward identity in the 
presence of a non-zero order parameter. This divergence occurs when we 
take a "twisted" quark mass to zero, and we conclude that quenched QCD 
with Wilson fermions is well-defined only with a non-vanishing twisted 
mass. We show that this peculiar behavior of the near-zero-mode den- 
sity is special to the quenched theory by demonstrating that this density 
vanishes in the unquenched theory outside the Aoki phase. We discuss 
the implications for domain-wall and overlap fermions constructed from 
a Wilson-Dirac kernel. We argue that both methods work outside the 
Aoki phase, but fail inside because of problems with locality and/or chiral 
symmetry, in both the quenched and unquenched theories. 



1. Introduction 



Wilson fermions play a prominent role in Lattice QCD pp. They (and their im- 
proved versions) are widely used in numerical calculations of hadronic quantities. In 
addition, they are at the heart of the construction of lattice Dirac operators with 
domain-wall |2l El or overlap jH Ej fermions. All these fermion methods based on the 
Wilson-Dirac operator preserve all of the flavor symmetry, but not ordinary chiral 
symmetry. For Wilson fermions, a tuning of the bare fermion masses is needed to 
restore chiral symmetry in the continuum limit, while domain- wall fermions (with 
infinite extent in the fifth dimension) and overlap fermions possess a modified version 
of chiral symmetry with essentially the same algebraic properties as the chiral sym- 
metry of the continuum theory [Hj. This lattice chiral symmetry reduces to that of 
the continuum theory in the continuum limit. 

In this paper, we will discuss the phase diagram in the gauge-couphng, quark-mass 
plane for two degenerate Wilson fermions, for both the quenched and unquenched 
theories. We will be interested in correlation functions constructed from the inverse 
of the two-flavor (hermitian) Wilson-Dirac operator, evaluated on an equilibrium 
ensemble of gauge-field configurations. The difference between the quenched and 
unquenched cases is that only in the unquenched case the fermion determinant (which 
is positive for two degenerate flavors) will be part of the Boltzmann weight used 
to generate the ensemble, whereas in the quenched case the Boltzmann weight is 
obtained from a local pure-gauge action. The quenched theory can be understood as 
a euclidean path integral with two physical quarks, and two "ghost" (bosonic) quarks 
with the same mass, whose role it is to cancel the physical-quark determinant jTj |H] . 

Long ago, a conjecture was made for this phase diagram jH], which is shown 
schematically in Fig. 1, where Qq is the bare coupling and tuq is the bare quark mass. 
We will be concerned with the region —8 < amo < because only in that region 
can the Wilson-Dirac operator have any zero or near-zero eigenvalues. Spontaneous 
symmetry breaking (SSB) occurs as a consequence of a non-zero density of such modes 
[TU] , and the interesting part of the phase diagram is thus confined to this region (the 
"super-critical" region). The usual continuum limit corresponds to the critical point 
at Qq = 0, arriQ = 0. Other continuum limits are obtained for any value of arrio by 
taking qq — > 0, but generically all quarks stay massive {i.e. have a mass of order of 
the cutoff). Only at —amo = 0, 2, 4, 6, 8 do massless quarks show up in the continuum 
limit, respectively 2, 8, 12, 8, 2 of them (in our two- flavor theory). The number of 
quarks is determined by the number of momenta (or "corners" ) in the Brillouin zone 
where the free Wilson-Dirac operator has a zero. 

According to this conjecture, in regions A and C no SSB takes place, while re- 
gion B, the Aoki phase, is defined by the existence of a pionic condensate, which 
breaks parity and flavor symmetry.^ Let us briefly review the evidence for the (un- 
quenched) Aoki phase diagram. First, a simple heuristic argument suggests that a 
pionic condensate must exist in some part of the phase diagram 0. Start with some 

"'^Note that all chiral symmetries are explicitly broken at non-zero lattice spacing, and thus play 
no role in determining the phase diagram. 
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Figure 1: A representation of the phase diagram proposed by Aoki for two-flavor 
QCD with Wilson fermions and standard plaquette action. The solid lines mark phase 
transitions (all believed to be continuous). Phase B is the Aoki phase, defined by the 
existence of a parity- and flavor-breaking pionic condensate. Its "trade mark" are the 
five "fingers" reaching to the critical points on the go = Q axis. Phases A and C have 
no condensate. The lightly shaded area marks the super- critical region where near- 
zero modes may occur. The lines tuq = and arriQ = —8, that define the boundaries 
of the super- critical region, appear to have no special dynamical significance: the A 
phases extend on both sides of them. The quenched phase diagram is discussed in the 
main text. The darker shaded area is roughly the area where domain-wall fermion 
simulations have been carried out. The dashed line represents a possible trajectory for 
taking the continuum limit with domain-wall or overlap fermions (see Sect. 7). 

mo > and examine what happens as rriQ is decreased. For a non-zero lattice spac- 
ing a {i.e. go > 0) there is no chiral symmetry. Instead of the continuum relation 

oc \mq\ (here m^r and rriq are respectively the pion and quark masses) we have 
that oc rriQ — m', where the fact that m' = m'{gQ) does not vanish is a conse- 
quence of the breaking of chiral symmetry on the lattice. For rriQ > m' the pions 
are massive, becoming massless at the critical line mo = m'{go). For mo < m', m^ 
would go negative, signalling the breaking of a symmetry. A pionic condensate forms 
in some direction in flavor space, and the line mo = m'{gQ) determines the location 
of a second-order phase transition. The corresponding pion becomes massive again 
for mo < m', while the other two pions become Goldstone bosons associated with 
the spontaneous breaking of the SU(2) flavor symmetry ("isospin") down to a U(l) 
symmetry. Since the condensate is pionic, it breaks parity symmetry as well. Micro- 
scopically, the condensate arises from near-zero modes of the Wilson-Dirac operator 
[TUj . and can thus only occur in the region —8 < amo < 0. 

This argument does not provide much information on the detailed form of the 
Aoki phase. Additional analytical evidence comes from several sources. The location 
of the critical points along the line = is obtained from weak-coupling perturbation 
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theory, which, however, gives no information on the existence of a condensate. The 
existence of a condensate was discussed in the context of the pion effective action 
in Refs. [TTl [T2j. Ref. |T2! showed that, z/ the Aoki phase extends all the way to 
(7o = 0, it does so as indicated by the "finger" structure^ in Fig. 1, with the width of 
these fingers proportional to {akQcoYi where Aqcd is the QCD scale. In the strong- 
coupling limit, the location of the two critical points was first found in Ref. pi], and 
the nature of the phase with broken symmetry was clarified in Ref. [Hj. 

The phase diagram was studied numerically in both quenched fSj and full PB] 
QCD. The numerical results provide evidence that the five critical points at (70 = 
are indeed continuously connected to the "main body" of the Aoki phase at large 
coupling. In the quenched theory, evidence for all ten critical lines at qq = 1 [jS = 6.0) 
was found in the last paper of Ref. ^Hj- In the present paper, in which we address 
the mechanism responsible for the existence of an Aoki phase, the detailed form of 
the phase diagram is not important. For the sake of argument, we will assume that 
the very plausible situation depicted in Fig. 1 describes the actual phase diagram. 
For a discussion expressing some doubts, see Ref. [T7j . 

A difficulty with the quenched simulations is that exceptional configurations 
were discarded. Exceptional configurations fHl OHj will play an important role below. 
Keeping this subtlety in mind, no light excitations were found in the quenched A and 
C phases, except close to the Aoki (B) phase boundary (and, in particular, close to 
the five critical points at go = 0). 

While this picture appears to be rather satisfactory, there exists some evidence 
that, in the quenched case, seems to disagree with the diagram of Fig. 1. Numerical 
studies of the near-zero modes of the (hermitian) Wilson-Dirac operator indicate 
that, for qq > 0, a non-zero density of near-zero modes always occurs in the quenched 
theory, anywhere in the super-critical region [201 ISl (namely, in regions B and C, and 
the super-critical part of region A). Through the Banks-Casher relation, this would 
imply that the pionic condensate vanishes nowhere in this region, and SSB takes place 
everywhere, thus contradicting the phase structure sketched in Fig. 1. Moreover, one 
would be inclined to expect Goldstone excitations everywhere in the super-critical 
region, in which case Fig. 1 would be completely wrong for the quenched theory; 
it would not even serve as a guide to the long-range physics. We note however, 
that this is in conflict with what is known from the numerical studies of Ref. [TH] 
mentioned above, as well as with an analytical study 22 , in which it is argued that 
the effective-field theory analysis of Ref. ^2] is also valid for the quenched case. 

Another clue comes from a study of a special class of zero modes by Berruto, 
Narayanan and Neuberger (which we will refer to as BNN hereafter) pH]. They 
showed that, for choices of mo in the super-critical region, the Wilson-Dirac operator 
has exact zero modes for very smooth gauge fields with one dislocation, which is 
contained in a small hypercube with a linear size of a few lattice spacings (later we 
give a more precise description of their results). These zero modes are exponentially 
localized, and their existence characterizes these "BNN" gauge-field configurations 
as exceptional configurations. They further argued that in large volume one would 

similar phase diagram is known to exist in the Gross-Neveu model 
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expect configurations with a "dilute gas" of sucli dislocations, because of their highly 
localized nature. This class of configurations could then contribute to a non- vanishing 
density of near-zero modes of the Wilson-Dirac operator for the quenched theory, 
where such configurations are not suppressed by the fermion determinant. This study 
thus lends analytical support to the numerical results reported in Refs. j^Hll^II- 

In this paper, we will propose a resolution to this apparently paradoxical situation. 
First, we prove that, if there is a non- vanishing density of exponentially localized 
near-zero modes, there exists a different mechanism for saturating the Ward identity 
involving the order parameter for SSB. The orientation of an order parameter is 
determined by a small external "magnetic field," here provided by a so-called twisted- 
mass term 121] • We argue that under the above circumstances, no Goldstone pole 
need appear, but instead that the two-point function of the would-be Goldstone-pion 
field diverges in the limit of vanishing twisted mass, even if the momentum does not 
vanish.^ This divergence already occurs in finite volume, and we will argue that it 
persists in the infinite- volume limit of the quenched theory. 

We then invoke a concept from condensed-matter physics in order to come up 
with a conjecture for the microscopic picture underlying the quenched phase diagram. 
The motivation is to understand what happens in the case that we have a dense, 
rather than dilute, gas of dislocations. For this purpose, we interpret the square of 
the hermitian Wilson-Dirac operator as the hamiltonian of a five-dimensional theory, 
and thus its (positive) eigenvalues as energy eigenvalues. The ensemble of gauge 
fields on which these eigenvalues are computed act as a random potential, and a non- 
zero measure subset of these fields can be viewed as random distributions of highly 
localized scattering centers. The eigenstates describe the possible states of a (four- 
plus-one-dimensional) "electron" in this background. Note that the fact that we are 
concerned with the quenched theory here makes the simplicity of this picture more 
compelling, because of the absence of any feedback of the fermions on the gauge fields. 

In general, there will exist a band of extended eigenstates above a certain energy, 
denoted A^, whereas eigenstates below this energy will be localized exponentially.'^ 
While attractive-potential scattering centers will tend to localize the electron as in 
a bound state of an isolated scattering center, a large-enough density of them will 
make it possible for the electron to travel throughout the lattice by tunneling, with 
the corresponding state becoming extended. The outcome will depend on the energy 
of the state and on the density and other properties of the scattering centers. The 
(average) localization range of the states will increase as the energy is increased 
towards A^, becoming infinite at A^, which in condensed-matter physics is referred to 
as the "mobility edge" ^ (for reviews see Refs. |2Hl 1211 EO] ) • 

The value of the mobility edge is a dynamical issue, which depends on the ensem- 
ble of gauge fields, and thus on the location in the phase diagram. We conjecture that 
the Aoki phase is precisely that region of the phase diagram in which the mobility 
edge is equal to zero. Outside the Aoki phase, it is larger than zero, and consequently, 

■^While completing this paper, we became aware of the fact that a similar observation has been 
made some time ago in condensed matter theory )25j . 

'''In a less detailed form, this physical picture has previously been considered in Ref. |2(ij . 
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all near-zero modes are exponentially localized. When one moves closer to the phase 
transition, the mobility edge comes down continuously, until it vanishes and the Aoki 
phase is entered. 

We complete our discussion of the physical picture by invoking our results on 
the two different mechanisms by which the Ward identity can be saturated in the 
presence of a non-zero condensate. When the mobility edge is zero, only extended 
near-zero modes contribute to the spectral density, and the Ward identity predicts 
the existence of Goldstone excitations that dominate the long-range physics of our 
theory; when the mobility edge is larger than zero, the condensate (and, hence, the 
order parameter) is produced by exponentially localized near-zero modes, the pion 
two-point function is diverging, and, we conjecture, no Goldstone excitations occur. 
Thus the region with Goldstone bosons coincides with the region where the mobility 
edge is equal to zero, and the existence of Goldstone bosons can be used as a definition 
of the Aoki phase in both unquenched and quenched QCD. 

This picture of what the quenched phase diagram looks like begs the question as 
to whether it is consistent with unquenched QCD. In the unquenched case, the Boltz- 
mann weight is modified by the fermion determinant, which tends to suppress the 
entropy of gauge-field configurations with near-zero modes. We show that the diver- 
gence in the pion two-point function cannot occur in the unquenched case. Therefore, 
indeed no near-zero-mode density can build up unless the resulting condensate is ac- 
companied by Goldstone bosons (in which case we are inside the Aoki phase), and 
exponentially localized near-zero modes do not play the same role in unquenched 
QCD as in quenched QCD. 

The scale of the typical localization length of the near-zero modes is set by the 
lattice spacing. The existence of a non-zero condensate outside the Aoki phase is thus 
a short-distance artifact. One expects no long-range physics (such as the existence 
of Goldstone excitations) as a consequence of localized near-zero modes, unless their 
average localization length becomes so large that they behave collectively. According 
to our conjecture, this happens precisely when the mobility edge comes down to zero, 
and, just as in unquenched QCD, only extended near-zero modes contribute to the 
spectral density. Therefore, an effective-lagrangian analysis in terms of the long-range 
effective degrees of freedom, as carried out in the unquenched case [12], should also 
make sense in the quenched case. Both inside and close to the boundaries of the Aoki 
phase, the effective lagrangian should provide a valid description of the long-range 
physics. Indeed, such an analysis has been carried out in the quenched theory 
and leads to conclusions very similar to those obtained in the unquenched case. 

Thus far, our results deal with the case of QCD with Wilson fermions. How- 
ever, there are important consequences for lattice QCD with domain-wall or over- 
lap fermions constructed from the Wilson-Dirac operator. In short, we claim that 
domain-wall and overlap fermions can only be used well outside the Aoki phase. For 
domain-wall fermions, the hermitian Wilson-Dirac operator is closely related to (the 
logarithm of) the transfer matrix that hops the fermions in the fifth dimension. A 
density of near-zero modes of the hermitian Wilson-Dirac operator implies the exis- 
tence of long-range correlations in the fifth direction. This threatens the decoupling 
of the left-handed and the right-handed quarks that live on opposite ends of the five- 
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dimensional world, and, hence, the restoration of chiral symmetry in the limit on an 
infinite fifth dimension. 

The crucial issue is the value of the mobility edge of the Wilson-Dirac operator. 
Well outside the Aoki phase the near- zero modes are all exponentially localized, with a 
scale set by the lattice spacing. As a result, their contribution to observables vanishes 
when the four-dimensional separation is set by a physical scale (without having to tune 
any parameter). As for extended modes of this operator, they cannot mediate long- 
range correlations in the fifth dimension, and thus they do not obstruct the recovery 
of chiral symmetry either. In contrast, inside the Aoki phase the near-zero modes are 
extended, and they mediate long-range correlations in all five directions. This strongly 
suggests that inside the Aoki phase the resulting four- dimensional effective theory 
either is completely non-local, or, at best, contains long-range degrees of freedom 
different from the desired ones. This physical picture is valid for any value of the 
lattice spacing in the fifth dimension as, and thus the same conclusion applies to 
overlap fermions which corresponds to the special case as 0. 

This paper is organized as follows. We define our theory, introduce the twisted 
quark mass, and derive the relevant Ward identity in Sect. 2, where we also briefiy re- 
view the relation between the condensate and the near- zero- mode density. In Sect. 3 
we derive a spectral representation for the pion two-point function. We observe that 
in the super-critical quenched theory this two-point function always diverges in finite 
volume, in the limit of a vanishing twisted mass. We thus provide an alternative 
mechanism for saturating the Ward identity in the presence of a non- vanishing con- 
densate. In Sect. 4, we start with defining localized near-zero modes more precisely. 
We show that a non-zero density of them will dominate the divergence of the pion 
two-point function in large volume, and, likely, will produce the divergence in the 
infinite-volume limit as well. We then give a qualitative but detailed discussion of the 
mobility edge, and its role in determining which way the Ward identity is satisfied. 
This leads to a comprehensive picture of the quenched phase diagram, and in partic- 
ular of the distinction of the super-critical regions inside and outside the Aoki phase. 
In Sect. 5 we discuss briefiy how the picture changes if we reintroduce the fermion 
determinant, i.e. if we unquench the theory. We then discuss the implications of our 
analysis for domain- wall and overlap fermions in Sect. 6, and finally summarize our 
conclusions in Sect. 7. In Appendix A we review Anderson's definition of localization 
as the (partial) absence of diffusion [27], and its relation to our definition in Sect. 4. 
Some technical details are relegated to Appendix B. 

2. Twisted-mass QCD, the Ward identity, and the Aoki con- 
densate 

The definition of the Wilson-Dirac operator is 




(2.1) 
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^^xy ^^xy ^ ^ ~l~ ^x—p,,yUy^ 

in which = (c?, i), where are the three Pauh matrices, and each entry is a 
2x2 matrix. The hermitian Wilson-Dirac operator is H{mo) = 75D(mo), with 
75 = diag (1,1,-1,-1). The theory is symmetric under the replacement amo 
— (8 + amo) because of the fact that for amo = —4 the Wilson-Dirac operator only 
contains nearest-neighbor couplings, thus allowing for a f/(l)e symmetry jl8j . 

We define HQ^mo) as the operator H{mo) with all U^fM = 1. This corresponds 
to the line 5^0 = in the phase diagram. The spectrum of Hlimo) covers a closed 
interval [(A^i„(mo))^ (A^„^(mo))2], with < A^i„(mo) < A[^„^(mo) < 00. Aj^^^ is 
determined by minimizing 



a^Hl{p] mo) = sin^(ap^) + ^(1 - cos {ap^)) + amo (2.2) 



over the Brillouin zone. Keeping three components of the momentum fixed, it is 
easy to see that a^HQ{p\mo) is linear in the cosine of the fourth component, and 
thus minimized when this cosine equals ±1. It follows that at a minimum all four 
components of the momenta have to be equal to or tt, and thus that a^ifg (p; mo) 
is the minimum over k of {2k + amo)^, in which k is the number of momentum 
components equal to vr. We thus find that A[J^j„ = min |mo — mgl, where am^ is one of 
the values 0, —2, —4, —6, and —8. We see that A[^j„ is generically of order 1/a, except 
when mo is close to one of these critical points. These critical points correspond to 
the tips of the Aoki "fingers" on the line (70 = in Fig. 1. Near these critical points 
the theory yields 1, 4, 6, 4, and 1 light quarks per Wilson fermion, respectively. All 
eigenmodes of H^lmo) are plane waves, and thus extend over the whole volume. Note 
that the spectrum of ifo(mo) is contained in two disjunct intervals, [—^maxy~'^rnin\ 
and [Aj^i„, X^,^^], separated by a gap (if A^^^ 0). 

Returning to the interacting theory, we consider for which values of mo the opera- 
tor H{mo), or equivalently D{mo), can have zero eigenvalues. Write D{mo) = A + iB, 
with A and B hermitian, and consider an eigenmode D"^ = [A + iB)'^ = A\E'. We 
then have that ^^{A - iB) = ^U*, and thus 2^t/i^ = (A + A*)^^^. It follows that 
A can only vanish if \E'M\E' vanishes. Since A = —{W + amo), this can only happen 
if — 8 < amo < 0, because the spectrum of W is confined to the interval [0,8]. The 
super-critical region {i.e. the region where zero modes may exist) is thus the region 
in which —8 < amo < 0, and we will restrict ourselves to that region for the rest of 
this paper. 

In the super-critical region "exceptional" configurations may, and do, occur. Ex- 
ceptional configurations are defined by the condition that D{mo) has an eigenmode 
\E'o with an exactly real eigenvalue Ao JHICSI- For such configurations we have that 
if (mo + \o)^ = 75-D(m.o -l- Xo)"^ = 0. Hence, a configuration is exceptional ijf H has 
an exact zero mode for some mo. BNN configurations (see Sect. 4) are a special kind 
of exceptional configurations. 
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Wc will be interested in a two-flavor theory constructed with this Wilson-Dirac 
operator, with a fermion lagrangian 



C = ijj{D{mo) - imi75T3)V' (2.3) 
= V {H{mo) - imiTa)^' , 

where ijj' = ijj^y^ is a field redefinition with jacobian one. The (eight-component) 
field t/j transforms in the fundamental representation of isospin SU(2). We added a 
symmetry-breaking term proportional to mi pointing in the T3 direction, where Tk is 
another set of Pauli matrices acting in isospin space. This symmetry-breaking term, 
which breaks both isospin and parity, and thus has the quantum numbers of one of 
the pions, allows us to probe the existence of an Aoki phase, in which this pion field 
develops a vacuum expectation value. We will assume the standard plaquette action 
for the gauge field, unless otherwise noted. Integrating over the fermion fields yields 
the fermionic partition function 

Zf = n(^n + ^ml)(A„ - im,) = ]l{Xl + ml) , (2.4) 

n n 

where the product is over the eigenvalues A„ of H{'mo), of which there is a finite num- 
ber on a finite-volume lattice. Note that Zp depends on the gauge-field configuration 
through the eigenvalues A„. 

The relevant Ward identity is obtained by performing a local fiavor transformation 

5+ip{x) — ia{x)T+ip{x) , 5+ip{x) — —ia{x)ip{x)T+ , (2-5) 

in which t± = (ri±ir2)/2. With n±{x) = iip{x)'y5T±ip{x) and rc3{x) = i'4'{x)j5T3ip{x), 
we find for any operator O that 

d; (J^ix) 0{y)) + 2m, Mx) 0{y)) = ^ {S^O{y)) , (2.6) 

(Jj 

where the backward lattice derivative is defined by d*f{x) = {f{x) — f{x — p,))/a and 
the corresponding vector current is 

•^m(^) = ^(V^(^)r+(7/. + l)C//.(^)V'(^ + /i)+V^(^ + /i)r+(7M-l)t^/I(^)V'(^)) ■ (2-7) 

While the notation (...) indicates an integration over both fermion and gauge fields, 
we note that the Ward identity is also valid if we integrate over the fermion fields 
only. Taking 0{y) = Ti-{y) and defining 

T{x,y) = (7r+(a;)7r_(y)) , (2.8) 
V,{x,y) = (j;(x)7r_(y)) , 

we arrive at the Ward identity 

a;r^(x, y) + 2m,T{x, y) = % (7r3(y)) . (2.9) 

a 
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Defining 



r(p) = ?7Ee^'^'"^^r(x,y), (2.10) 

xy 



and similarly r^(p), the Fourier-transformed Ward identity is 

i ^(1 _ e-^-r>,^V,{p) + 2mif(p) = (773) . (2.11) 

Tlie pionic condensate 

(vr3) = (aVV^)E(^3(y)) (2.12) 

y 

can be calculated by first calculating the expectation value of vr3(?/) in a fixed gauge 
field (denoted by (. . .)u) from Zp, 



a' E {My))u = Z^'ji-Zp = 2 E , (2.13) 



mi 

and then averaging this over the gauge field, to obtain 



(7r3)=2/rfAp(A)^^^, (2.14) 

where p(A) is the spectral density defined below. In the limit mi — > 0, we obtain 

(7r3) = 27rp(0) , (2.15) 

the Banks-Casher relation ^] for the case at hand. 

We define the spectral density from the cumulative density 

a^(a) = (e , (2.16) 

\a„<a / 

where \I'„(x) is the eigenfunction associated with the eigenvalue A„, 

H{mo)^n = Xn^n, (2.17) 

and with normalization a'^J2x\'^n{x)\'^ = 1. Because of the gauge-field average in 
Eq. ()2.16|) the right-hand side is actually independent of x. The spectral density p(A) 

p{X) =dAf/dX = V'^ ^n)^ ■ (2-18) 

All these results are valid in finite- volume QCD, in both the quenched and unquenched 
cases. They remain true if the infinite-volume limit is taken, but of course this limit 
may or may not commute with the limit mi — > 0. It is well known that in unquenched 
QCD these limits do not commute in a phase with SSB; the order parameter (7r3) in 
the mi = theory vanishes in finite volume, but does not vanish if the thermodynamic 
limit is taken before the limit mi 0. We will return to this point in Sect. 5. 
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For any given mo we will denote by A^j„ i^max) minimum (maximum) eigen- 
value of H'^imo) over the gauge-field space. These values are finite, because of the 
fact that H{mo) connects a given site only to a finite number of neighboring sites, 
and because the gauge group is compact; hence if (mo) is uniformly bounded. The 
cumulative density ^f{\), defined by the spectrum of H^ttiq), is monotonously non- 
decreasing on the interval [—Xmax, ^max]- Using the results of Ref. jSHj, we will argue 
below that A^m = for any > amo > —8 and that, in the quenched theory, p(A) > 
for any —Xmax < A < Xmax- We will first consider a finite volume (Sect. 3) and then 
the infinite-volume hmit (Sect. 4). Outside the super-critical region, i.e. for amo > 
or amo < ~8) one has Xmin > 0, and p(A) = for —Xmin < X < Xmin- 

3. Spectral representation of the pion two-point function 

We are now ready to establish one of the key results. Using the spectral repre- 
sentation for the quark propagators, it is straightforward to derive that, in a finite 
volume. 



V{x,y) = f VUB{U)Y.'^i{x)^,{x) \ ^\{y)^M . V • (3-1) 
J , Ai. -|- ^'m^ Ar, — ^m^ 



The two propagators on the right-hand side correspond to the "up" and the "down" 
quarks respectively. Each sum runs over the eigenstates of the (single-fiavor) her- 
mitian Wilson-Dirac operator H{mo). In full QCD the Boltzmann weight is ;B = 
Zpe^pi—Sg), where Sg is the gauge action, and the fermionic partition function Zp 
is given in Eq. ()2.4|) . In quenched QCD, B = exp{—Sg). 

For a generic gauge field, the eigenvalue spectrum of H{mo) is non-degenerate: 
if An = Afc then "^nix) = \E'fc(x). Gauge field configurations for which this is not true 
form a subset of the configuration space with measure zero, because either if(mo) 
would have to have a symmetry in this fixed gauge-field background, or there is an 
accidental degeneracy. Both types of degeneracy disappear under a generic small 
deformation of the gauge field. 

In a finite volume on the lattice, the sums over n and k in Eq. ()3.1|) are finite, and 
the integration over gauge fields is absolutely convergent because we integrate over 
a compact group. This implies that we may interchange the sums and the integral, 
and consider separately the contribution for each n and k to r{x,y), which is 



For n ^ k (and thus generically A„ 7^ A^), expression ()3.2j) is finite, even in the limit 
nil ^ 0. This is because the denominators A„ + irrii and A^ — irrii become singular 
on different subspaces (with codimension one) of the gauge-field space, and the inii 
terms in the denominators provide an ie prescription for how to integrate around the 
poles at A„ = or Afc = 0. Therefore, in the limit mi only terms with n = k can 
make a non-zero contribution to the product mir{x,y) : 



kn 




(3.2) 



hm mir(x,i/) = lim / VUBiU)Y.\^M\^\^^{y)\^ 

ni— >0 mi^O J 



X^ + m' 



2 • 

■1 



(3.3) 
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Of course, the fact that we keep the volume finite is a key element of this argument. 
For the Fourier transform of r{x,y), Eq. (j2.im) . this result translates into 

lim miT^p) 

mi^O 

Hnip) 

Eq. ()3.3p has a dramatic consequence for quenched QCD. We may introduce the 
density correlation function 

n{x,y;X) = /l?Wi3(W)^|^„(x)n*„(7/)|2<5(A-A„), (3.5) 
with which Eq. ()3.3|) can be written as 

dXTZ(x,y; X)— ^ = 7r7^(x, 0) . (3.6) 

X'' + mf 

It is clear that lZ{x,y; X) is independent of mi in the quenched theory because B{U) 
is. Also, if we choose mo anywhere inside the super-critical region, the existence of 
BNN-like zero modes^ implies that there is no spectral gap around A = 0, and, in 
particular, that exact zero modes occur on a non-empty subspace of codimension 
one. As a result, 7l{x,y; A) > for any —Xmax < X < Xmax, because, apart from the 
delta function, the integrand in Eq. ()3.5j) is strictly positive. Hence, we find a 1/mi 
divergence in T{x,y), 

7r7^(x,?/;0) 

^{x,y) = + 0(1). (3.7) 

mi 

Moreover, if 7l{x, y; 0) > 0, also p(0) > and, from Eq. ()2.15|) . (tts) > for the same 
reasons, and we find that we can have SSB in a finite volume in the quenched theory.® 
Eq. (j3.7p shows that the quenched theory in the super-critical region is singular in 
the limit mi — * 0, and it provides an alternative mechanism for saturating the Ward 
identity Eq. ()2.9p . With this alternative mechanism, no Goldstone excitations need 
to appear in the quenched theory, even in the presence of SSB. 

The fact that these results are so simple is related to the fact that, in the finite- 
volume quenched theory, the existence of both the condensate and the 1/mi diver- 
gence in the pion two-point function are kinematical effects, which depend only on 
the super-criticality of ttiq and the strict positivity of the (mi-independent) quenched 
measure. The magnitude of (tts) > and 7l{x,y] 0) depend on the dynamics, and, as 
it will turn out, on whether we are inside or outside the Aoki phase. 

To conclude this section, we note that one can repeat the entire analysis assuming 
a different SU(2)-orientation of the flavor- and parity-symmetry breaking term in the 
lagrangian (j2.3j) . The orientations of both the pionic condensate and the divergent 
term in the pion two-point function follow the orientation of this symmetry breaking 
term, just as in the case of conventional SSB. In the next section, we will address the 

^For a more detailed account of BNN's work, see Sect. 4. 
^For the unquenched theory, see Sect. 5. 



iu„!;ii/OT)i:^. (3.4) 

mi^o V J „ Xi + mi 



V„ + ml 

a^5:|vl/„(x)|V^P^ 
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dynamics of the quenched theory, and in the section after that the differences with 
the unquenched theory. 

4. Localization and spontaneous symmetry breaking in the 
quenched theory 

In a given finite volume, all eigenmodes of if(mo) are localized for the trivial 
reason that their support is compact. In particular, all near-zero modes contribute 
to the 1/mi divergence we found in the previous section. (Here we define a near- 
zero eigenvalue as an eigenvalue with an absolute value of order mi, for a given 
(small) mi.) However, as we will now demonstrate, quantitatively the finite- volume 
divergence comes from exponentially localized near-zero modes. We postpone the 
discussion of what happens in the infinite-volume limit to the second part of this 
section. 

4-1 Localized modes and the divergence of the pion two-point function 

On a finite lattice of spacing a and (large) linear size L (in all four directions), a 
normalized eigenstate "^nix) is exponentially localized, provided there exist a positive 
constant Ci = 0(1), a lattice site (which depends on \l/„), and some / with a < 
I <^ L such that 

|vl/„(x)p< ^exp(^-^^j . (4.1) 

The localization range (or localization length) In is the minimal value of / for which the 
bound 1)4.11) holds. For definiteness, the distance \x — yl"^ is defined by the minimum 
of J2ij, — + rif^Ll^ over all integers n^. Our intention is to establish a lower 
bound on T{p) in terms of the density of localized near- zero modes and, for this 
purpose, we take Ci to be independent of / and L. Also, the restriction to /„ < / <^ L 
means that we consider only localized modes whose support (of roughly order Z^), and 
size ~ inside this support, are basically independent of the volume. This 

definition may exclude some modes which would be exponentially localized according 
to some other reasonable definition, but the class we are considering here will be 
sufficient for our purpose. 

We begin with establishing bounds on Tinip) iii Eq. (|3.4p : 

l-c,llp'<\nn{p)?<l. (4.2) 

Here C2 oc ci is another numerical constant which we will define below. The upper 
bound is trivial, while the lower bound just expresses the fact that one cannot resolve 
the structure of a localized mode with localization range In using momenta p <C 1/Zn. 
The lower bound is established by noting that (with the "center" of the localized 
eigenstate of Eq. ()4.1|) ) 

0<1-Re(e*^'<7^„(p)) = 2a^^sin2(p(x-x°)/2)|^„(x)p (4.3) 

X 
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_ C{ln) ,22^^ ,2 2 

The inequalities ()4.1|) and | sin(a)| < a were used. The dimensionless quantity 

C(0-4i:^=^-'1._, (4.4) 

is a continuous function of / and finite in the hmit I ^ oo. The fact that C{1) is finite 
for / — i> oo is a consequence of choosing the bound on with the prefactor 

Ci/l^, with Ci independent of /. This is important for deriving a non-trivial lower 
bound on T{p). It thus follows that 

C2 =m.axi>a{C{l)} (4.5) 

is a finite numerical constant of the same order as Ci. Since \Re {e'P< nn{p))\ < 
\T-Cn{p)\ the inequality ()4.2|1 follows. 

We now use this bound to derive a lower bound on T{p) in Eq. ()H.4j) . for mi 0. 
Choosing some fixed /, we split the eigenstates of if (mo) into two classes: those which 
are exponentially localized according to our definition of Eq. 1)4.11) with localization 
range /„ < /, and the rest. Keeping only the former, we obtain 

mif (p) > ^ / VUB{U) Y: K^'f'^f + 0{m^) . (4.6) 

Since in this sum we only keep terms for which In ^ this may be expressed more 
simply as 

mif(p) > npi{0) (l - C2 ly) + 0(mi) . (4.7) 

Here pi{X) is the density of exponentially localized eigenstates (according to our def- 
inition ()4.ip ) with eigenvalue A and localization range less than or equal to /. 

This result is of key importance. It establishes that if there is a non-zero density 
of near- zero modes which are exponentially localized within a certain range I, the 
pion two-point function diverges as mi in any finite volume, for all momenta p 
below (roughly) the inverse localization range. As we will argue in Sect 4.2 below, 
the analysis of Ref. allows us to establish the same result in the infinite-volume 
limit, provided the mobility edge does not vanish (and mo is in the super-critical 
region). This provides an alternative mechanism for saturating the Ward identity, 
Eq. ()2.11|) . in the presence of a non- vanishing pion condensate (tts). It also shows 
that the quenched theory in infinite volume is not well defined for mi = 0, at least 
as long as the mobility edge does not vanish. 

A comment is in order on the range of momenta for which our result is valid, 
because of the fact that in a finite volume the momenta are discrete, and thus cannot 
be chosen arbitrarily small. First, we observe that for a localized near-zero mode, 
the scale of the typical localization range I is set by the inverse mobility edge, A~^. 
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(This will be explained in more detail in the following subsection.) As long as we are 
not close to the Aoki phase, this is of order one in lattice units. With the smallest 
lattice momenta being of order we may choose the volume large enough that 
C2{1/LY ^ 1, so that non-zero momenta indeed exist for which the right-hand side 
of Eq. ()4.7|) is positive. 

4-2 Spontaneous symmetry breaking without Goldstone excitations 

We are now at a point where we wish to collect all the available results, and 
use them to construct a conjecture for the quenched phase diagram. The goal is 
to combine the numerical E] as well as the analytical |i2^ evidence that basi- 
cally everywhere in the super-critical region the density of near-zero modes does not 
vanish (already in finite volume) with the inequality for the pion two-point function, 
Eq. ()4.7p . In a sense, we will be proposing the simplest possible picture of what deter- 
mines the phase diagram, while requiring consistency with the Ward identity ()2.11|) . 
While we already summarized our conjecture for the phase diagram in the introduc- 
tion, we present a more detailed discussion in this subsection. 

A key element is the observation made by BNN in Ref. [22 5 that in the super- 
critical region there will be a density of near- zero modes. We will therefore start with 
elaborating on this observation. As already alluded to in the introduction, it will lead 
us to borrow the concept of the "mobility edge" from condensed-matter physics as 
the crucial ingredient in determining the quenched phase diagram. 

Following BNN, we begin with an infinite- volume configuration containing a single 
dislocation. We pick a lattice gauge field which is equal to the classical vacuum 
{Ux^l = I) everywhere, except for the links of an rt^ hypercube containing the origin. 
These links can take any value. The square of the Wilson-Dirac operator can be 
written as 

H\m^)=Hl{m^) + V , (4.8) 

where ifo(^o) is the free operator. The potential V is supported on an {n + 2)^ 
hypercube. The "hamiltonian" H'^{mo) is a discretized, semi-positive Schrodinger 
operator with a finite-range potential, with a continuum threshold at (A[^j„)^ (c/. 
Sect. 2). All eigenstates with (A^j„)^ < E < (A[J^^^)^ are scattering states, while 
those with < E < (A^j„)^, if they exist, are bound states of if^(mo) with "binding 
energy" 

-E,^E-iXl,J<0. (4.9) 

Of course, bound-state wave functions decrease exponentially away from the hyper- 
cube, with a decay rate determined by E^-^. In particular, a zero mode has E = 0, 
and thus a binding energy Ef, = (Aj^j„)^. 

For single- dislocation configurations, BNN proved that the problem of finding a 
zero mode in large (or infinite) volume can be reduced to a corresponding eigenvalue 
problem on the small hypercube. Solving the latter problem numerically BNN then 
found that, for — 1 > arriQ > —7 (roughly), the links of a hypercube as small as 2^ can 
always be chosen such that iJ(mo) has an eigenvalue equal (or numerically extremely 
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close) to zero. This means that, for that range of ttiq, zero modes exist for any volume 
down to 2^. Moreover, these zero modes are exponentially localized (as long as mo is 
not equal to one of its critical values). We note that BNN configurations are examples 
of exceptional configurations.'' 

It is not implausible that, by allowing for larger dislocations, exponentially- 
localized zero modes would be found throughout virtually the entire super-critical 
interval > amo > —8. One does expect that, to get a zero eigenmode closer to a 
boundary point of the super-critical interval, the gauge field would have to differ from 
the classical vacuum over bigger hypercubes. In this paper, we will assume that this 
is the case, i.e. that for any arriQ strictly inside the interval (—8,0), BNN-like zero 
modes will be found. 

Before continuing, let us comment on the nature of typical zero modes, and the 
configurations supporting them. An exact zero mode of if^(mo) is also a zero mode of 
H{mo). As one varies the lattice gauge field, a zero eigenvalue generically corresponds 
to a zero-level crossing, because with Wilson fermions there is no chiral symmetry 
that might protect a zero eigenvalue from moving away from zero. In other words, if 
W(r) = {Ux^{t)} is a family of gauge-field configurations that depend smoothly on 
r , and if there is an eigenvalue such that A„(W(0)) = then, typically, A„(W(r)) 
changes sign at r = 0. This is true for (almost) all exceptional configurations which 
support exact zero modes. It implies that any finite-volume configuration supporting 
an (exponentially localized) exact zero mode has some open neighborhood, in the 
gauge-field space, with a spectrum of near-zero modes. Also, in a finite volume, the 
subspace of configurations that support an exact zero mode has codimension one. 

Our next step is to consider configurations containing a very small density of 
BNN-like dislocations. The dislocations are small, surrounded by the classical vac- 
uum, and (still, on average) far apart. This is a "controlled" form of disorder. The 
origin of disorder is two-fold: the positions of the dislocations are chosen at random 
(one can speak about a dilute gas of dislocations); also the links that define each dis- 
location are chosen at random. A unique limiting value lim|2.|^oo V{x) of the potential 
in Eq. ()4.8|) no longer exists, because dislocations may be found arbitrarily far from 
the origin. 

As long as the gas of dislocations is dilute, we can still identify all dislocations 
separately from the surrounding vacuum, simply by inspection. So, let us focus on 
one dislocation, and assume that in isolation {i.e. with no other dislocations any- 
where), it has an exact zero mode. Since the other dislocations are far apart, they 
will have a negligible effect on this zero mode. This is easily seen by invoking a 
variational argument. Using the eigenfunction of the zero mode of the isolated dislo- 
cation as a trial wave function, the expectation value of if^(mo) on this trial state will 
receive contributions only from the other distant dislocations. If the mean distance 
between dislocations is R, the expectation value of H'^{mo) will be of order e~'^'^ where 
7 ~ A[^j„i?. Therefore H'^{mo) must have an eigenstate with an exponentially small 
eigenvalue E ~ e~'^'^ which closely resembles the original zero mode. We may now 

^Even though we expect BNN configurations to amount to only a tiny subset of all exceptional 
configurations. 
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vary the gauge field at the dislocation such that the eigenvalue under consideration 
varies as well. As long as R is large enough (compared to 1/A^j„), e~'^'^ is small 
enough that this gauge-field variation will vary the eigenvalue over an interval that 
includes zero. 

The situation changes qualitatively if the density of dislocations is large, and also 
will be different for configurations generated in a typical Monte-Carlo simulation. 
In order to describe this situation, it will be useful to introduce the mobility edge 
hypothesis |2Zj. The mobility edge is well-defined only in infinite volume, and we will 
thus assume the volume to be infinite for the following discussion. The hypothesis 
states that, when disorder is introduced in an ordered system, the conduction-band 
structure of the ordered system is replaced by a number of alternating energy intervals, 
each containing either (exponentially) localized or extended eigenstates. In particular, 
no energy exists for which there are both localized and extended eigenstates.^ The 
energy separating an interval with localized modes from one with extended modes is 
referred to mobility edge. 

Intuitively, a mobility edge arises as follows. Let us first return to a situation with 
"controlled" disorder consisting of a dilute gas of BNN dislocations. As already men- 
tioned, for a very small density of dislocations, all bound states with Ei, = 0(l/a^) 
(see Eq. ()4.9j) ) will be exponentially localized near a single dislocation. As the den- 
sity of dislocations is increased there is an increased probability that an "electron," 
situated initially in a bound state of a given dislocation, will tunnel into a near-by 
dislocation. When the probability to tunnel a certain distance times the average num- 
ber of available dislocations within this distance has become 0(1) for some given E, 
the electron will be able to travel infinitely far by tunneling from dislocation to dis- 
location, regardless of the details of the dislocation that produced the original bound 
state. In other words, the eigenstate with eigenvalue E becomes extended. 

For a single dislocation, any eigenstate with E < (Aj^j„)^ decays exponentially. 
Increasing E means a slower decay rate and, hence, an increased tunneling probability. 
This means that eigenstates with higher E become extended first when the density 
of dislocations is increased. If is the mobility edge, the range < E < will 
consist of exponentially localized eigenstates, while the range E > consists of 
extended ones. For a very small density of dislocations the mobility edge will be close 
to (Aj^j„)^, while for larger densities it will move further away from this value. If, for 
instance, we allow only dislocations generating an attractive potential in Eq. (|4.8|) . 
the mobility edge will clearly be lower than (Aj^j^)^. 

We note that there may exist more than one mobility edge at a given location 
in the phase diagram. For instance, if the hamiltonian is bounded from above, one 
expects another mobility edge to occur, above which other localized eigenstates exist. 
In our case, we expect a second mobility edge like this, because the eigenvalue spec- 
trum of the free hamiltonian, HqIttlq) in Eq. ()2.2j) . also has a maximum eigenvalue 
iKnaxY- thi^ paper, we will be concerned with the lowest-lying mobility edge, and 
in particular, the question of whether it vanishes or not. 

Now consider the hypothetical situation that, for a certain gauge- field configura- 

^We are not aware of any proof of this — widely used — assertion. 



17 



tion, both extended and localized eigenstates of the hamiltonian H^Iitlq) occur at the 
same energy E. If we then calculate the eigenstates on a typical small fluctuation of 
this gauge field, the extended and localized states will mix, and all new eigenstates 
will be extended. The only possible exception is when all eigenstates of the original 
gauge field were localized. A "typical" configuration thus has either only extended 
or only localized eigenstates at a given E. Since typical configurations (in this sense) 
determine the properties of an ensemble, it follows that the value of a mobility edge 
separating localized and extended eigenstates is associated with an equilibrium en- 
semble, and thus with a point (mo, go) in the phase diagram. Another way of arriving 
at the same conclusion is to note that in the thermodynamic limit, one equilibrium 
configuration suffices to determine the equilibrium properties of the theory, and thus 
is typical. For our purposes, it will be useful to characterize each point (amo,(y'o) 
in the phase diagram by the value of the lowest mobility edge of if^(mo). (In the 
rest of this paper, we will refer to this mobility edge as "the" mobility edge.) The 
mobility-edge hypothesis will be at the heart of our conjecture for the phase diagram 
of quenched QCD. 

Before we get to this conjecture, let us briefly comment on the average localization 
range as a function of E. For energies at which only localized states exist, one expects 
to be able to define an average localization range I = 1{E). This average localization 
range should then grow if the energy E approaches the mobility edge, and diverge at 
and beyond the edge, where only extended eigenstates exist. 

Based on BNN's results we expect that anywhere inside the super-critical region 
there will be a non-zero density of near-zero modes. If the mobility edge is zero this 
expectation is fulfilled by assumption. If the mobility edge is larger than zero, then, 
for any > aruQ > —8 and go > 0, there is a finite probability per unit volume 
that any given Monte-Carlo configuration will contain a ball with (a large) radius R, 
inside which the configuration will resemble the classical vacuum, apart from a single 
BNN dislocation well inside that ball. Using the same variational argument as in the 
dilute-gas case, the exact BNN zero mode associated with the dislocation implies the 
existence of an eigenstate with eigenvalue e~'^'^ -C A^, well below the mobility edge. 
This eigenstate decays exponentially both inside and outside of the classical-vacuum 
ball. By invoking small deformations of this configuration, the existence of a density 
of exponentially localized near-zero modes follows.^ 

We now formulate our main conjecture about the phase diagram of quenched 
QCD with two Wilson fermions. Near the line go = 0, we expect the mobility edge to 
be very close to (A^j„)^. In particular, with (A[J^j„)^, it will vanish at the critical "end" 
points go = 0, amo = 0, —2, —4, —6, —8. As go is increased and the typical gauge-field 
configuration in an equilibrium ensemble becomes less smooth, the mobility edge will 
move away from (A^^^)^, and, in parts of the phase diagram, it may vanish. Since 

zero-level crossing in the spectrum of (mo) changes the index of the overlap operator. This 
index was advocated as a definition of the topological charge on the lattice However, standard 

arguments imply that the topological charge density scales like 1/VV in large volume, and, thus, 
vanishes for V ^ oo (see e.g. Ref. (^). Therefore, topological considerations do not explain a 
non-zero density of near-zero modes in the infinite volume limit. This situation reflects the inherent 
ambiguity in trying to identify small dislocations as instantons or as anti-instantons. 
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it vanishes at the critical end points on the go = axis, it is reasonable to expect 
that the mobility edge will also vanish in the vicinity of those points for qq > 0, 
thus opening up the finger-like structure of region B in Fig. 1. In a region where the 
mobility edge vanishes, the density of near- zero modes will be due to extended near- 
zero modes, and one expects long-range behavior, in particular Goldstone excitations. 
We thus conjecture that, in quenched QCD, the Aoki phase — which is defined as the 
phase in which Goldstone excitations exist — coincides with the region of the phase 
diagram where the mobility edge vanishes. The near-zero mode density, and thus the 
value of {tts), is due to extended eigenstates. Outside the Aoki phase (but inside the 
super-critical region), the mobility edge is larger than zero, and the density is due to 
exponentially localized near-zero modes. Because of this non-vanishing density, the 
mechanism derived in the previous subsection will kick in, and the pion two-point 
function will diverge for mi — > for small-enough momenta. 

Extended states should not contribute to the 1/mi divergence in the infinite- 
volume limit. The number of extended states grows like the volume V, whereas 
the contribution to Eq. ()3.3p of each extended state drops like This yields a 

contribution of order 1/V, which vanishes for V — ^ oo. Therefore, we expect no 1/mi 
divergence of the pion two-point function inside the Aoki phase; the 1 / mi divergence 
characterizes region C (and the super-critical part of region A) of the phase diagram. 
It follows that, while (tt^) is not a useful order parameter for detecting the Aoki phase, 
the residue of a 1/mi divergence in T{p) (for small enough momentum) is a useful 
order parameter in this sense. There is even a local order parameter associated with 
this divergence: 

$^[x) = lim mir{x,x) = lim mi (7?+ (x) 7r_ (x 

mi—*0 mi— >0 

This follows from ^(x) = (l/V) J2p^{p), the bound ()4.7|) . and the positivity of T{p) 
for all p. 

We believe that essentially the same conjecture holds for unquenched QCD. The 
only difference is that in the unquenched case, the spectral density p(0) (and thus the 
order parameter) vanishes outside the Aoki phase, because of the suppression of the 
localized near- zero modes by the fermion determinant (c/. Sect. 5). Therefore, also 
the 1 / mi divergence in the pion two-point function will not occur in the unquenched 
case. These phenomena are quenched artifacts, but the existence and role of the 
mobility edge are not. To summarize, the qualitative features of the phase diagram 
in Fig. 1 are valid in both quenched and unquenched QCD. The mobility edge Ac is 
zero in region B, and non-zero outside of this region, where exponentially localized 
modes with < A^ exist. In regions C and in the super-critical part of regions A 
the spectrum of localized modes extends down to zero, with the only difference that 
p(A) > for A^ > in the quenched theory, and for A^ > in the unquenched theory. 

Clearly, our picture of the phase diagram proposed here is a conjecture, and we 
have no proof that it is correct. We emphasize that it appears to be the simplest 
possible way one can understand the phase diagram, given the available analytical 
and numerical evidence. We will end this section with a more quantitative, but still 
heuristic argument as to why Goldstone excitations only occur in the phase with 
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(4.10) 



vanishing mobility edge. 

We will assume that outside the Aoki phase, the near-zero-mode density p(0) is 
entirely due to localized zero modes captured by our bound, Eq. (|4.ip . Moreover, we 
will assume that, for practical purposes, a maximum localization length l^ax exists 
for the near-zero modes. Hence^° Pimaxi^) — sjid, using both the lower and 

upper bounds in Eq. ()4.2|) . 

7rp(0)(l - < hm mif(p) < 7rp(0) . (4.11) 

mi— >0 

Writing r^(j9) = ip^T^ijP')^ it then follows from the Ward identity, Eq. (j2.11|) that 

S(0) = 27rp(0)0(O, (4.12) 

leading to the conclusion that T^{p) does not have a massless pole. 

The most drastic assumption we make here is that there is a finite Imax- This 
is rather unlikely, even if it is reasonable to expect that the probability to find a 
near-zero mode with a very large localization range at some point well outside the 
Aoki phase is very small. In order to argue that also in this case no Goldstone poles 
occur outside the Aoki phase, clearly one needs a better estimate of |H„(p)| than that 
provided by the bounds given in Eq. ()4.2j) . While it is a very hard problem to come 
up with better estimate, we will use an ansatz for a better estimate to show how an 
improved argument might work. Let us assume that |7i„(p)| can be estimated by 

I^"(P)I~T-^' (4-13) 
1 + pH^ 

This is not unreasonable for momenta p ^ 1/a. For very small momenta, one expects 
only the long-distance features of the near-zero modes, i.e. the exponential tail in 
Eq. (j4.ip . to determine the Fourier transform of |\l/„(x) and thus YKnip) I would have 
to look something like this. The precise form is not important; the crucial properties 
of this ansatz are that it is bounded for arbitrary /„, and that it is consistent with 
the fact that |7i„(0)| = 1, as follows from the fact that \E'n(a;) is normalized to one. 

If p(0) is completely due to localized near-zero modes, we may write p(0) = 
/(f d/p;(0) with p'i{X) = dpi{X)/dl, and Eq. Q as 

27r r°° r 1 
lim2mir(p) ~ --/ dl VUB{U)Y.K\Ml-Q (4.14) 

POO I 

Combining this with Eqs. ()2.1H) and ()2.15p we obtain the estimate 

^ip')-2nl dim^^^^,, (4.15) 

showing again that T^{p) has no massless pole. We note that the integral over / 
should be convergent because dpi{0)/dl should be strongly suppressed for large /. 

^"Recall that pi{X) is the density of locaUzed eigenstates with locaUzation range less than I. 
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This completes our discussion of the quenched phase diagram. In the remainder 
of this paper, we will expand on the differences between the quenched and unquenched 
cases, and explore the consequences for domain-wall and overlap fermions. 

5. Unquenched QCD 

In this section we only review well-known facts about unquenched QCD. The 
only reason we include this brief review is to point out that none of the surprising 
results we obtained in the previous section for quenched QCD are in conflict with the 
standard lore in the unquenched theory. 

5. 1 Vanishing of the condensate in finite volume and absence of localized 
near-zero modes 

We start with a review of the proof that the finite-volume condensate is zero 
in unquenched QCD. Summing the Ward identity ()2.9j) over space-time, the total- 
derivative term drops out, and we obtain 

2mia^^(7r+(x)7r_(y)) = (7r3) . (5.1) 

X 

In a finite volume, the correlation function of the product of any (finite) number of 
fermion (and link) operators is bounded. The left-hand side of Eq. ()5.1|) is therefore 
the product of mi with a bounded function, hence it vanishes for mi — > 0. It follows 
that the condensate vanishes. This is the familiar result that there is no spontaneous 
symmetry breaking in a finite volume. 

In order to prove boundedness, consider first the unnormalized expectation value 
of any observable made of a product of (link and) fermion fields. Because of the 
Berezin integration rules for Grassmann variables, the Wick contraction of the fermion 
fields leads to a function which is analytic in the parameters of the fermion action 
mo and mi. Integrating also over the (compact) link variables of the gauge field 
leaves an unnormalized expectation value which again is an analytic function of these 
variables, including for mi = 0. Similarly the partition function Z is strictly positive 
for a two-flavor theory with Wilson fermions (for qq > 0), and again is bounded 
for a bounded range of values of tuq and mi which includes mi = 0. Hence normalized 
expectation values are bounded as well. In contrast, in quenched QCD the fermion 
determinant is missing, and a fermionic observable may diverge if the Wilson-Dirac 
operator has near- zero modes. This is what happens in the case of the pion two-point 
function. 

It follows, through the Banks-Casher relation, that p(0) = 0. Since < pi{0) < 
p(0), one has pi{0) = for any finite volume. While of course the infinite volume is 
subtle (and does not commute with the limit mi — 0), we believe that for any finite /, 
P/(0) will remain zero in the infinite- volume limit, because pi{0) reflects short-distance 
physics. Also, if the inequality ()4.7|) remains valid in the infinite-volume limit on the 
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one hand, and unquenched QCD does not have a 1/mi divergence on the other hand, 
clearly pi{0) must remain zero in that limit. The fact that a 1/mi divergence indeed 
cannot occur in unquenched QCD is explained in the next subsection. 



5.2 The Goldstone theorem 

Here we re-derive the Goldstone theorem in the euclidean path-integral context 
in infinite volume (for a review see Ref. |S2])- For simplicity, we do this in the 
continuum; nothing relevant changes on the lattice. As in the previous subsection 
our aim is to show that the (Fourier-transformed) pion two-point function T{p) is 
bounded in a neighborhood of mi = but now for p 7^ (and without assuming 
that the condensate necessarily vanishes). We may then again conclude that, for 
p ^ 0, 2mir{p) for mi — > 0. Taking this limit, the Ward identity ()2.11|) yields 
(dropping terms of order ap) 

ip^f^(p) = (vrg) , p^O. (5.2) 

If the condensate is non-zero {i.e. if SSB takes place), this implies the existence of a 
massless Goldstone pole in T^{p). 

To prove the boundedness of T{p) for p 7^ 0, we invoke the Kallen-Lehmann 
representation. A one-particle contribution to T{p) must be of the form {p'^+M'^) 
where M is the particle's mass and F is a form factor which is non-zero on shell 
(p2 = -M^). For any non-zero euclidean momentum, 

mi mi , , , 

Py^O. (5.3) 



p2 + M2 - p2 

The right-hand side of inequality ()5.3|) provides a bound for the contribution of any 
one-particle excitation to 2mir(p), implying that it vanishes in the limit mi — > 0, 
for p ^ 0. This is true also if the particle's mass vanishes in the same limit (as in 
the case of a Goldstone boson). The contributions of multi-particle states will be 
less infra-red singular, and should also vanish. In summary, in unquenched QCD, 
2mir(p) vanishes in the thermodynamical limit, for p 7^ 0. 

The quenched theory is not unitary, and the Kallen-Lehmann representation for 
the pion two-point function invoked above is not valid. Apparently, nothing stops the 
pion two-point function from developing a 1/mi divergence. 



6. Implications for domain- wall and overlap fermions 

The most common constructions of lattice Dirac operators with domain-wall or 
overlap fermions employ the Wilson-Dirac operator discussed in this paper (or im- 
proved versions thereof) ClS Sb kernel. 

Domain-wall fermions are five-dimensional Wilson fermions, in which only hop- 
ping terms in the four physical directions couple to four-dimensional link variables, 
which themselves are independent of the fifth coordinate [21 E] . In the most common 



22 



version, the fifth dimension is restricted to a finite interval of length a^Ng, where 
and Ng are respectively the lattice spacing and number of sites in the fifth direction. 
Free boundary conditions are employed on either side (in the limit of a vanishing phys- 
ical quark mass). In the limit Ng oo, massless four-dimensional fermions appear 
which are bound to the two boundaries. If the left-handed component of this massless 
fermion is bound to one boundary, the right-handed component is bound to the other 
one. For qq ^ 0, precisely one such massless fermion appears (per five-dimensional 
fermion field) if the domain-wall height M = —mo is chosen such that < a'M < 2 
where a' = max{a, 05} [33,. This means that niQ has to be super-critical. 
The overlap-Dirac operator is defined as j3] 



with iJ(mo) the hermitian Wilson-Dirac operator. Notice that 7I = 1. For this 
operator to describe one massless flavor in the continuum limit, one needs to choose 
> aniQ > —2. 

Because of the fact that if(mo) plays a crucial role in the construction of domain- 
wall/overlap fermions, it is natural to expect that the phase diagram of QCD with 
Wilson fermions has important implications for properties of domain- wall/overlap 
fermions. The most important dynamical issue is whether the mobility edge is (close 
to) zero, or not. In this section, we will argue that in order to retain locality in lattice 
QCD with overlap fermions, the parameters of the lattice theory must be chosen well 
outside the Aoki phase. For domain-wall fermions the situation is equivalent: only 
well outside the Aoki phase will chiral symmetry be restored exponentially fast with 
increasing A^^- In addition, only in that case will the four-dimensional effective theory 
(which is described by a generalized overlap operator, the details of which depend on 
05) be local in the limit Ng ^ 00. 

In the quenched case, we may say that lattice QCD with domain- wall/over lap 
fermions is inside or outside the Aoki phase if the "underlying" theory with the 
Wilson-Dirac operator H{mo) is. One may think of the phase diagram as that of 
a theory with Nf quenched domain- wall/overlap flavors and two quenched Wilson 
flavors. Of course, this changes nothing in the correlation functions of the domain- 
wall/overlap quarks. What it means is that, also in the domain- wall/overlap case, we 
deflne the quenched Aoki phase by the existence of a Goldstone pole in the appropri- 
ate correlation function constructed from the inverse of the Wilson-Dirac operator 
if (mo) (c/. Eq. ()2.8p ). using the same value of mo as in the kernel of the domain- 
wall/overlap operator, and the same Boltzmann weight. According to our conjecture, 
this corresponds to the region in parameter space in which the mobility edge of if (mo) 
vanishes. 

In the unquenched case, the only thing that changes is the Boltzmann weight 
used to generate the ensemble of gauge-fleld conflgurations on which the domain-wall 
or overlap operator is computed. For any given ensemble, the mobility edge of 

^^Note that the fermion determinant used for generating the gauge- field configurations may corre- 
spond to any type of lattice fermion. If the sea quarks and the valence quarks are not coming from 
the same fermion action, one is dealing with a partially quenched theory |34|. 
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the Wilson-Dirac kernel H{mQ) should have a well-defined value. We may thus still 
define the Aoki phase as that region in parameter space in which the mobility edge 
of Himo) vanishes. Conversely, if the mobility edge does not vanish, we will say that 
we are outside the Aoki phase for the unquenched theory. 

Let us start with locality of the overlap operator. The overlap operator Dov{x, y) 
cannot have a finite range j^, and thus is not strictly local for a finite lattice spacing. 
However, Hernandez, Jansen and Liischer proved that the overlap is local in the sense 
that \Dov{x,y)\ decays exponentially with the distance \x — y\, provided the gauge 
field obeys an admissibility condition [36 . The effect of this condition is to secure 
an 0(l/a) gap in the spectrum of Hijuo) (except very close to the critical values of 
mo). They furthermore generalized this result to the case that if (mo) has an isolated 
zero mode inside an otherwise empty spectral gap by showing that such a mode is 
necessarily exponentially localized. The rate of the exponential decay of the overlap 
operator in these cases is of order one in lattice units. 

The problem is that, for realistic simulations, it is impractical to impose an ad- 
missibility condition on the gauge fields. If instead one uses one of the commonly used 
local gauge actions to generate an equilibrium ensemble, H{mo) will have (localized) 
eigenmodes with very small eigenvalues, and these eigenvalues will not be isolated in 
the sense of Ref. jHEI- In fact, since the number of localized modes grows in propor- 
tion to the volume, in the infinite- volume limit the eigenvalues of the localized modes 
will form a dense set, and no eigenvalue will be isolated. This is most clear in the 
quenched case, where p(0) is always non-zero in the super-critical region. 

At this point we invoke our physical picture of the phase diagram. For any qq > 0, 
the band edge A[^j„ (of the free hamiltonian Ho{mo)) gets replaced by the mobility 
edge Ac. We hypothesize that, as long as Ac > 0, the conclusion of Ref. jSHI still holds: 
if all near-zero modes are exponentially localized with a finite average localization 
length I, then 75 (and thus Dov) decays exponentially. As we argue below, the decay 
rate is in principle governed by the smaller of Ac and l/I. We expect this situation to 
apply inside phase C in Fig. 1, well outside the Aoki phase. In fact, if we assume, like 
at the end of Sect. 4.2, that the near-zero modes have a finite maximal localization 
length Imax outside the Aoki phase, our hypothesis is simple to prove. 

First, however, let us consider the corresponding situation with domain- wall 
fermions. The relevant Ward identity [Sj for the study of the finite- A^^ chiral symmetry 
breaking with domain-wall fermions is, for x ^ y, 

d;{A^,i^) Jl{y)) = 2m,{j^{x) Jl{y)) + 2{jt,{x) Jl{y)) . (6.2) 

Here is the domain-wall partially-conserved axial current and J5 is the correspond- 
ing pseudo-scalar quark density {a,b label fiavor generators); rriq is the quark mass. 
Jgq is a pseudo-scalar density located midway between the boundaries of the five- 
dimensional bulk. The term containing this density represents the chiral symmetry 
breaking at finite Ng other than the expected breaking coming from an explicit quark 
mass. For chiral symmetry to be restored, it should vanish in the Ng ^ 00 limit for 
non-singlet axial currents. (In the case of the axial U(l) current this term gives rise 
to the anomaly in the continuum limit.) 



24 



In order to study chiral symmetry restoration, it is useful to consider the transfer 
matrix T{M, 05) which hops in the fifth direction, from one four- dimensional "time" 
slice to the next OH]- For every eigenmode \E'„ of T(M, 05) with (positive^^) eigen- 
value Un, we let Qn = min{ujn,uj~^}. Hence, < g„ < 1. In the second-quantized 
transfer matrix, the replacement of Un by g„ reflects normal ordering. For the free 
theory the transfer matrix has a gap: g„ < go for all eigenmodes, where go < 1- 
We refer to go as the band edge of the free transfer matrix. (As an example, in 
the special case where aM = a^/a = 1, one has go = 1/2.) In the interacting the- 
ory, we define the mobility edge of the transfer matrix as gc = max{g„}, where the 
maximum is taken over the extended modes only. Below, we will also speak of the 
"hamiltonian" H{—M, 05) = — log {T{M, a^))/a^. This hamiltonian has a band edge 
— Iog(go)/a5 in the free theory, and a mobility edge = — log(gc)/a5 in the inter- 
acting theory. (In the limit 05 — »• one recovers the familiar Wilson-Dirac operator: 
H{-M,0) = H{-M), and reduces to A^.) 

Since the transfer matrix has a gap in the free theory, the four-dimensional mass- 
less fermions at the boundaries are exponentially bound to their respective boundaries. 
If we take A^^ ^ cxo at fixed M and 05, we find an exponentially decreasing overlap 
between the (fifth-coordinate) wave functions of the light-quark modes tied to the 
boundaries and the "midway" pseudo-scalar density Jg^. Hence chiral symmetry gets 
restored exponentially in Ng. With gauge fields obeying an admissibility condition, 
the situation remains the same, because the transfer matrix still has a gap. Also, 
with commonly used local gauge actions, a similar behavior has been demonstrated 
in weak-coupling perturbation theory ^Tf IHHj. 

Non-perturbatively, one can prove that, for any 05, a zero mode of the Wilson 
kernel is also an eigenstate with eigenvalue one of the transfer matrix: T(M, a5)\E'n = 
^„ if and only if i7(-M)^„ = H- This implies that H{-M) and H{-M,a^) 
have a similar spectrum of near-zero modes. A near-zero mode of H{—M, a^) causes 
long-range correlations in the s direction. This means that in realistic simulations 
with commonly used gauge actions, configurations will occur for which the last term 
in Eq. ()6.2|) may be large, thus "threatening" the chiral symmetry of domain-wall 
fermions even in the large- A^, limit. This danger is again brought out most clearly 
in the quenched case, where a non-zero density of near-zero modes always occurs for 
the above specified range of M. 

However, we hypothesize that, in analogy with the overlap case, there is a funda- 
mental difference between the effect of exponentially localized near-zero modes and 
that of extended ones. If all near-zero modes of H{—M, 05) are exponentially localized 
with an average localization length 7, we expect their contribution to the symmetry- 
breaking term in Eq. (j6.2|) to decay exponentially as a function of the four- dimensional 
separation \x — y\. If the scale of I is set by the (four-dimensional) lattice spacing, 
this decay will resemble a contribution from excitations with mass of the order of the 
cutoff. When the lattice spacing is small enough, even at finite A^, this contribution 

^^For a certain range of M and as the transfer matrix may have real but negative eigenvalues; in 
this case the entire analysis can be carried out in terms of the transfer-matrix squared (provided Ns 
is a multiple of four) . 
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will vanish for large |a; — ?/| relative to the other two terms in Eq. ()6.2|) . whose long- 
distance behavior is determined by a physical mass. All other contributions to the 
symmetry-breaking term will not be suppressed with the four-dimensional separation, 
but they will vanish exponentially with Ng. 

Similar conclusions apply to the locality of the effective four-dimensional lattice 
Dirac operator in the Ng —>■ oo limit |^]. This is a generalized overlap operator 
constructed by making the replacement H{—M) —>■ H[—M, 05) in Eq. (jfj.lj) . Because 
of the similar zero-mode structure of H[—M) and H{—M, 05), we expect the effective 
Dirac operator to be local (in the exponential sense) well outside the Aoki phase, as 
in the case of the overlap discussed above. 

For both domain-wall and overlap fermions, these arguments break down if the 
near-zero modes are extended, or even, if l/I becomes of the same order as the nu- 
merical values of the physical masses in a particular simulation. This latter situation 
is expected close to the Aoki phase, because I will increase going towards the phase 
transition, and become infinite at the phase transition. The domain-wall formal- 
ism gives a clear indication on what may go wrong inside the Aoki phase. When 
the mobility edge of the Wilson kernel H{—M) (or, more generally, of H{—M,a5)) 
is zero, there are massless excitations everywhere inside the five-dimensional bulk. 
Normally, the contribution of the massive bulk modes is canceled by Pauli-Villars 
(pseudo-fermion) fields. But there is absolutely no guarantee that that cancellation 
will persist when the bulk fermion and pseudo-fermion modes are massless: this situ- 
ation corresponds to a different phase of the theory. The limit Ng ^ 00 now involves 
infinitely many unphysical, light four- dimensional fields (arising from both the five- 
dimensional domain- wall and the pseudo-fermion fields) , and we have every reason to 
worry that the limiting theory is not what we want it to be. This is true for any a^, 
and thus also in the overlap limit 05 — 0. 

We now present a more detailed argument as to why the overlap operator should 
be local well outside the Aoki phase. We will assume, as we did at the end of Sect. 4.2, 
that at a given point well outside the Aoki phase there exists a finite maximum 
localization length for the near- zero modes of H{mo). In fact, to avoid technical 
complications, it is convenient to make a slightly stronger assumption, namely, that 
if Ac > is the mobility edge, then all eigenstates with —Ac/2 < A < Ac/2 have 
a bounded localization length / < Imax < 00. As before, this assumption is not 
unreasonable because the average localization length diverges only for |A| y Ac, and 
the probability for arbitrarily large localization lengths is expected to vanish rapidly 
(and uniformly) for any |A| < Ac/2. By simply replacing if (mo) by if(— M, 05), the 
below argument may also be applied to the Ng ^ 00 limit of domain- wall fermions. 

We begin by noting that the overlap operator Dov is local if 75 is local. We split 

75 = 75<+75>, (6.3) 

where 7^ is the projection of 75 onto the subspace spanned by eigenstates with eigen- 
value |A| < Ac/2. On the basis of eigenmodes of H{mQ), 7^ can be represented 

^•^For small 05 one can also use Borici's kernel jlOl l5D{—M){2 — a5D{—M))^^ that gives rise to 
the same generalized overlap operator as H{—M, 05). 
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as 

75<(x,y) = a^ E ^n(x)^vl/t(^). (6.4) 

|An|<A,/2 

Assuming that all modes with |A| < Ac/2 satisfy inequality ()4.1|1 one has 

|75<(x,y)| < C,a' E ^ l \'^-'^n\^ \y-<\ \ (g_5^ 

|A„|<A,/2 V / 

Performing the ensemble average as in previous sections (see in particular Eqs. (j.'Kij) 
and dmH)) we find 

/ V /"CO /'Ac/2 

(l75<(x,y)|)<ciaM rf/ / d\p[{X)K{\x-y\ll), (6.6) 

where 

ir(|x - y|//) = l-^ I d^xo e-(l--ol+ls/-o|)/(20 ^ e-|x-,|/(2/) _ ^g^^) 

The last approximate equality holds for |x— y| ^ /, where the integral is dominated by 
localized modes supported inside a tube of radius I around the straight line connecting 
X and y. We now invoke the assumption that In < Imax for all —Ac/2 < A„ < Ac/2. 
The "worst-case scenario" is that the /-integral in Eq. ()6.6p will be dominated by 
/ ^ hnax- In this case we find the bound 

(|75<(x,i/)|) < exp(-|a; - y\/{2l^ax)) , (6.8) 

where constants and power corrections have been ignored. Hence 7^ is local. 

The operator 75* lives in the orthogonal subspace, whose projection operator is 
(x, y) = 5.j,^y -V<{x,y). Here 

V4x,y)=a' Yl ^n(x)vl>t (y) = (75<)^ (6.9) 

|A„|<Ac/2 

is the projector on the eigenstates with |A| < Ac/2. The bound ()6.5j] evidently applies 
to V<^{x,y) as well. Therefore 'P< is local, and, hence, also P> is local. Proceeding 
exactly as in Ref. [32] we write 

^> = H{mo)\H{mo)\-'V:.. (6.10) 

(Strictly speaking, the right-hand side is defined via its mode expansion.) The ex- 
pansion of (the eigenvalues of) \H{mQ)\~^ in terms of Legendre polynomials may now 
be invoked, and the locality of 75* follows. 

Let us discuss some implications of this result. In principle one can envisage 
two extreme situations. If Imax <S A~^ then the non-locality of 7^ and V> may be 
neglected. The localization range of the overlap will be governed by the mobility 
edge Ac, and will be of order A~^. In other words, in the absence of an admissibility 
constraint, A^ plays the role of the lower bound on the spectrum of if^(mo). The 
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opposite extreme is that l^ax ^ where we expect the locahzation range of the 
overlap to coincide with Imax- 

In our argument, we did not attempt to maintain the same level of rigor as in 
Ref. |36J. The advantage of our more heuristic argument is that it deals with the more 
realistic case of a density of localized near-zero modes where the methods of Ref. 
are inapplicable. Our analysis is only semi-realistic because, in principle, there will be 
a very tiny, but non-zero probability to encounter an exponentially localized near-zero 
eigenmode with an arbitrarily large localization length (in a large enough volume). In 
practice, however, we expect that our assumption of the existence of a finite, not too 
large Imax will be valid for most simulations. Of course, the most important question 
in any given simulation is how the scale of the typical localization length as well as 
the size of the mobility edge compare quantitatively to the scale of the physics one is 
trying to compute. 

We now turn to do main- wall fermions with finite A^^- A useful measure of chi- 
ral symmetry breaking is the so-called residual mass rrires IS], which is essentially 
the ratio of the two correlators on the right-hand side of Eq. ()6.2|) . Denoting the 
(euclidean) time coordinate by r we define 



where for convenience we have switched to flavor- changing densities. (In terms of the 
"isospin" symmetry of the two flavors occurring in these densities, the ± superscripts 
correspond to the operators r± of Eq. ()2.5j) .) 

The study of rrires consists of two steps. First, using the transfer matrix formalism 
of Ref. we derive expressions for the correlators in a fixed gauge-field background. 
For the numerator in Eq. ()6.11|) this step is outlined in Appendix B. Next we have to 
carry out the integration over the gauge field. In general this is quite complicated, 
and so we will restrict ourselves to two relatively simple cases. Observing that r 
is a measure of the four- dimensional separation in Eq. (jfj.llll . we will discuss the r 
dependence for fixed A^s ^ 1, and the Ns dependence for fixed |r| ^ a. As before, 
well inside the C phase we expect the mobility edge of H{—M, as) to be a quantity 
of order one in lattice units. For r = 0(a) and A^^ 3> 1, rrires will be dominated by 
the exponentially localized near-zero modes, whereas for A^^ = 0(1) and |r| ^ a it 
will be dominated by the extended states close to the mobility edge. 

The crucial observation is that, as explained earlier, the quark and anti-quark 
operators of the pseudo-scalar density J^(0, 0) live on the two boundaries of the fifth 
dimension, whereas J^q{x,T) is located midway between the two boundaries. In a 
diagrammatic language, the fermion (anti-fermion) at the boundary must be con- 
nected by a fermion line to the "midway" anti-fermion (fermion). In Eq. (jB.lll . these 
two fermion lines come from the one-particle sector of the second-quantized transfer 
matrix that acts Ns/2 times on the fermion (or anti-fermion) at the boundary. The 
contribution of each mode to a fermion line involves a factor of \E'ra(0, 0) q^"^'^ ^n(^? ''")? 
or its hermitian conjugate. 

For |t| > a the denominator in Eq. (jb.lip will contain short-distance contributions 




(6.11) 
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which are of no interest to us. In this case it is simpler to consider the numerator 
alone. Assuming some fixed A^^ ^ 1) the contribution of all extended modes has died 
out, because for them, Qn"^"^ '^Ic^"^'^ 1- We are left with the contribution of the 
exponentially localized modes of H{—M, a^), of which the near-zero modes dominate 
for A^s ^ 1. Following closely the overlap case, Eq. ()4.1|) allows us to put a bound on 
^'„(0, 0) "^n^x, t). a similar bound applies to the contribution of the second fermion 
line. The product of the two bounds gives rise to the same exponential factor as on 
the right-hand side of Eq. except that the factor of 1/2 inside the exponent is 

missing. As before we now assume the existence of a maximal localization length for 
the near-zero modes, arriving at 

{j+{x, r) J^{6, 0)) ~ exp(-|r|/U,) , AT, » 1 . (6.12) 

X 

The evident analogy with Eq. ()6.8|) is not surprising. 

Finally, we consider the A^^ dependence for |r| ^ a. Returning to rrires, now all 
the exponentially-localized modes may be neglected. We expect that rrires will be 
dominated by the extended modes of the transfer matrix with eigenvalues close to 
the mobility edge. The anticipated result is 

rrires ~ Qc" ) |r| > a. (6.13) 

This refiects a cancellation of the r dependence between the numerator and the de- 
nominator in Eq. ()6.1H) . for large r separations (for numerical evidence supporting 
this, see Ref. 01]). We are unable to derive this result analytically. We will give a 
heuristic argument, based on an analogy with perturbation theory, which supports 
this result. 

In perturbation theory, this result is obtained as follows jHH]- When the quark 
mass is zero, the one-loop domain-wall fermion propagator has, near each boundary, 
a factorizable form. Assuming for definiteness that the right-handed quark resides 
near the s = boundary, the propagator in the vicinity of this boundary is 

G{x, y; s, s') ^ x{s) Pr G(x, y) x\s ) . (6.14) 

Here Pr = |(1 + 75), x,y are the usual four-dimensional coordinates, and s,s' are 
coordinates in the fifth dimension. The separation |x — ?/| is assumed to be large com- 
pared to the lattice scale (but small enough for weak-coupling perturbation theory 
to be applicable). G{x,y) is the effective quark propagator whose tree-level Fourier 
transform is 1/^ for small p, and which, at higher orders, contains the familiar loga- 
rithmic self-energy corrections of a massless quark. The s-coordinate wave function 
x{s) carries no Dirac indices. 

For an optimally-chosen M as a function of go, the (tadpole- improved) tree- level 
wave function Xoi^) is completely confined to the boundary layer s = 0. For s > 1, 
the wave function arises from a short-distance, one- loop quantum effect. Explicitly, 

X{s) = Xo{s) + gl Xquantum{s) whcrC Xo(s) = (5s,0 ^Xid Xquan tum{s) ~ Up tO a pOWCr 

correction. This form of Xquantum{.s) gives rise to Eq. ()6.13|) . with go- (Since go is 
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the band edge of the free transfer matrix/^ the shape of the leading-order Xquantum{s) 
depends on the free domain-wall fermion action only. In higher orders one expects 
X,quantum (s) ~ iqof{go)Y where f{go) = 1 + 0{g^), and where f{go) depends on the 
gauge action too.) 

Restoring the momentum dependence, the s-coordinate wave function xi^jP) is 
universal in the sense that corrections to x{s) = x{s] 0) vanish like (ap)^. The physical 
reason is that the propagation in the fifth direction is dominated by a small region 
in the Brillouin zone surrounding the critical momenta that saturate the band edge 
go- These momenta are 0(1/ a). Fluctuations of the gauge field allow low- momentum 
modes to couple to these lattice-scale modes, with an amplitude that can be naturally 
expanded in powers of (ap)^. Eq. ()6.14|) is obtained by keeping only the leading term 
in this expansion. For all other terms, the 1/^ singularity of the (free) propagator is 
wiped out, and so their effect is negligible for large space-time separations. 

Non-perturbatively, we may envisage a similar factorization of the domain-wall 
fermion propagator for |r| ^ a (in a fixed gauge). Again, the extended states close 
to the mobility edge of the transfer matrix will mediate the s-propagation. Again 
these extended states will be controlled by the lattice scale, and it is reasonable that 
their coupling to low-momentum external states will have a universal value. Turning 
to the correlators of the gauge-invariant pseudo-scalar densities, a factorizable form 
of the (fixed-gauge) domain- wall fermion propagator, as in Eq. ()6.14j) . implies that 
the common r dependence of the numerator and the denominator of Eq. ()6.11|) is 
exp(— |r|m,r)- The result is Eq. (I6.13|l .^^ 

7. Conclusions 

Let us summarize our conclusions and make some additional comments. We start 
with what we learned about the quenched phase diagram. 

In finite volume, we argued that the Aoki condensate is non-zero everywhere in 
the super-critical region of the phase diagram, and that the pion two-point function 
always has a 1/mi divergence. This 1/mi divergence arises mainly from exponentially 
localized near-zero modes. If the restricted spectral density pi{0) is non-zero for some 
localization length /, then the momentum-space pion two-point function T{p) exhibits 
the 1/mi divergence for all momenta up to p^/^ ~ 1. Because of this divergence, 
clearly the finite- volume quenched theory is only well defined with a non- vanishing 
twisted mass mi. 

Extending this to infinite volume, we argued moreover that if all near- zero modes 
are localized, the condensate is approximately equal to limmi^o 2mir(p), and the 
difference vanishes with p^. This implies that there are no Goldstone excitations. 

Adopting the mobility-edge hypothesis we arrive at the following physical pic- 
ture. In the super-critical quenched theory the pionic condensate is always non-zero. 

^''in Ref. PE], the band edge of the free transfer matrix was denoted qi (and not go)- 

^^For the DBW2 gauge action at quenched ~ 2 GeV it was found that Qc ~ 0.6 021- This is 

very close to the value go = 0.5 found in one- loop perturbation theory [HH]- Therefore, in this case 

both higher-order effects and non-perturbative effects are small. 
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Goldstone poles exist, however, only in part of the super-critical region which, by 
definition, is the quenched Aoki phase. Inside the Aoki phase the mobility edge van- 
ishes, all the near-zero modes are extended, and the pion two-point function has no 
1/mi divergence, because the contribution of extended modes to this divergence goes 
to zero in the infinite- volume limit. Outside of the Aoki phase the mobility edge is 
larger than zero, and all the near-zero modes are exponentially localized. There are 
no Goldstone poles, and the Ward identity ()2.9j) is saturated instead by the 1/mi 
divergence of the pion two-point function. 

We can, in fact, completely characterize the quenched phase diagram using two 
local order parameters. For the SU(2)-orientation of the twisted-mass term in the ac- 
tion ()2.3p . one order parameter is the usual condensate {tc-^) = 27rp(0). The other is ^, 
defined in Eq. ()4.1()j) . which measures the size of the 1/mi divergence in the pion two- 
point function. In the quenched theory (tts) 7^ 0, and SU(2) is spontaneously broken, 
in the entire super-critical region. Hence the lines mo = and arrio = —8 are phase 
boundaries in the quenched theory. Inside the super-critical region we have a phase 
with Goldstone bosons (the Aoki phase), and phases with no Goldstone bosons (the 
C phases, and the super-critical parts of the A regions, which form separate phases 
in the quenched theory). The order parameter for having no Goldstone excitations is 

e>o. 

In continuum, infinite- volume, quenched chiral perturbation theory, it is known 
that the usual condensate (^''4'') develops a log(m), or with 6 small, divergence in 
the limit m — > 0113 • In the continuum limit of the lattice theory chiral symmetry 
is restored, and the Aoki condensate can be rotated back to the usual condensate 
(see for example Ref. [12] )• Turning this argument around, one expects a log(mi) 
divergence of the Aoki condensate inside the Aoki phase in the continuum limit, 
with 0{ami) corrections due to a non-zero lattice spacing. As explained above, we 
expect no 1/mi divergence of the pion two-point function inside the Aoki phase in 
the infinite- volume limit. 

In unquenched QCD we do not have any new results. We did verify that, only in 
unquenched QCD, one can prove that the pion two-point function T{p) is bounded. 
In finite volume this is always true, while in infinite volume this is true provided 
p 0. Thus, our new results for the quenched theory do not contradict any of the 
well-known facts for the unquenched case. 

The main differences between unquenched and quenched QCD are as follows. 
Inside the Aoki phase, the differences are accounted for by the chiral lagrangian, 
and the quenched theory is known to have enhanced chiral logarithms. Outside the 
Aoki phase, both unquenched and quenched QCD have exponentially localized near- 
zero modes; in the infinite- volume limit, zero is not an isolated eigenvalue of any 
gauge-field configuration drawn according to the Boltzmann weight. However, in the 
quenched theory p(0) 7^ 0, whereas in the full theory p(0) = 0. Heuristically, this 
can be understood as follows. Let us compare two infinite-volume gauge-field con- 
figurations, which differ only inside a small region whose radius is 0(1) in lattice 
units. Assume that, inside this small region, each configuration supports an expo- 
nentially localized mode, with possibly different eigenvalues. In the quenched theory, 
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these two configurations can have similar Boltzmann weights; but in the unquenched 
two-flavor theory, each locahzed mode will contribute to the Boltzmann weight of the 
corresponding conflguration (and, hence, to the spectral density) a factor of + mf 
{cf. Eq. ()2.4j) ). Therefore the configuration with the smaller eigenvalue will be sup- 
pressed. In the limit mi — 0, a configuration supporting an exponentially localized 
zero mode is completely suppressed in full QCD. In this argument, it is important 
that the difference between the two configurations is confined to a small region. If all 
the near-zero modes are extended, a local change in the configuration will have only 
an infinitesimal effect on any one of the eigenvalues, and the entropy may be large 
enough to overcome the suppression factor due to the fermion determinant. Thus, 
the extended modes may still build up a condensate, despite the suppression by the 
fermion determinant. 

On the basis of our analysis, we argued that the domain-wall and overlap formu- 
lations work only inside the C phase (s).^^ Inside those phases the mobility edge is 
non-zero, and all the near-zero modes of the Wilson kernel are exponentially localized. 
The overlap should be local, and for domain-wall fermions, chiral symmetry should be 
recovered exponentially with the size of the fifth dimension. Also, the effective four- 
dimensional overlap operator emerging in the Ng —>■ oo limit (which in the 05 — > 
limit coincides with the standard overlap constructed with a Wilson-fermion kernel) 
should be local. Inside the Aoki phase one encounters infinitely many light unphysi- 
cal modes, which contribute to the logarithm of the partition function with opposite 
signs. The overlap operator corresponds to the special case 05 0. It is hard to 
imagine how either formulation could remain valid under these circumstances.^^ 

It is worth commenting on the role of the gauge action in this respect. In par- 
ticular, let us discuss how precisely an admissibility condition jSHI would change the 
picture. An admissibility condition means that all the plaquette variables are con- 
strained to be closer to one than some < e <^ 1. The spectrum of if^(mo) then 
has a lower bound (A^^^)^ — 6^ (again A^^^ refers to the free hamiltonian ifo("^o))) 
where 6 is determined in terms of e. This prevents the existence of near-zero modes 
when A[^j„ > 6, but not when A^j„ < 6. The super-critical region will now be de- 
fined as that region in which H{mo) has zero modes when restricted to gauge fields 
satisfying the admissibility condition. This new super-critical region consists of five 
(disconnected) vertical strips in the phase diagram, each of which is located near one 
of the five critical value of rriQ. Inside each of those super-critical strips there should 
still be an Aoki phase. 

It is clear that, in principle, an admissibility condition guarantees the locality of 
the (ordinary or generalized) overlap, and the exponential recovery with A^^ of the 
domain-wall's chiral symmetry, for mg not close to a critical value. However, imposing 
an admissibility condition in numerical simulations is prohibitively expensive. Very 

^^One quark field per one lattice- fermion field is obtained by taking the continuum limit inside 
the C phase that borders on the interval —2 < amo < on the go — axis (where M — —mo is the 
domain-wall height, and provided < 1 33 ). A possible trajectory for taking the continuum limit 
is indicated by the dashed line in Fig. 1. 

-•^^For other work pointing at difficulties at strong coupling, see Ref. 01]. For related work on the 
phase structure of overlap fermions with a small hopping parameter, see Ref. |45| . 
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similar results can however be achieved by modifying the lattice gauge and/or fermion 
action, such that the density of the localized near- zero modes of if (mo) decreases for 
a fixed value of the lattice spacing. The density of near-zero modes was studied 
directly and indirectly in numerous publications; for recent reviews see Ref. |46j . 
In an approximation where all the exponentially localized modes are neglected, the 
mobility edge turns into a gap, which can be studied in perturbation theory |HH] 
(see also Ref. (SEl)- Non-perturbatively, numerical results with the Iwasaki and 
DBW2 gauge actions show a dramatic depletion of near-zero modes. This likely 
corresponds to larger values for the mobility edge at fixed lattice spacing. If we would 
replace the vertical axis in Fig. 1 by the lattice spacing (in physical units) this would 
amount to a recess of the Aoki-phase boundaries towards larger a, i. e. to enlarged C 
phases. 

We have studied the PCAC (Partially- Conserved Axial Current) relation for 
domain- wall fermions, Eq. (j6.2j) . The residual mass defined in Eq. (j6.1H) directly 
provides information on the rate of chiral symmetry restoration as a function of Ng. 
We argued that, by monitoring the anomalous term in the domain- wall PCAC re- 
lation, one can determine two crucial features of the "hamiltonian" H{—M, 05) that 
controls the propagation in the s-direction: the mobility edge (equivalently qc), cf. 
Eq. (jnHSl), and (in effect) the average localization length I of the near-zero modes, cf. 
Eq. ()6.12|) . These results may be combined as 

g-l^lA 

where c is a constant. (Here we have traded the maximal localization length Imax 
of the near- zero modes, assumed in Sect. 6, by an average localization length I.) 
The function /(r) contains the r dependence of the domain-wall's pion two-point 
function in the denominator in Eq. (jfj.llj) .^^ This result should be valid provided 
Ns'> 1 and/or |r| > a. 

The residual mass has been extensively studied in domain-wall fermion simula- 
tions which according to our terminology correspond to the range \t\ ^ a.^^ As a 
function of A^^, the results for the Wilson gauge action at quenched P = 6.0 and op- 
timally chosen M jll] are characterized by a rapid initial drop of rrires- However, for 
larger values of Ns the fall-off slows down, and eventually rrires settles at a non-zero 
value. Therefore, the first term on the right-hand side of Eq. ()7.H) is non-negligible 
in these simulations. According to our physical picture, this suggests that both the 
density of near-zero modes and their average localization length may be relatively 
large; hence, in this case, one is very close to the Aoki phase. In contrast, the 

-'^^As explained in sect. 6, for |r| > a Eq. 1)7.1(1 is not directly amenable to numerical tests because 
of the unknown short-distance effects in the pion two-point function; it may be advantageous to 
study to numerator of Eq. (|6.11(l directly, cf. Eq. (|6.12() . 

""^^The simulations often show a nice plateau mresiT, Ng) ~ mres{Ns) for a wide range of r. 

^°We believe that the first exploratory finite-temperature domain- wall simulations 0H| were in fact 
carried out inside the Aoki phase, which would explain the poor chiral symmetry observed even at 
relatively large Ng. This is why, in Fig. 1, we have drawn the region where domain- wall simulations 
have been tried such that it partially overlaps with the Aoki phase. 
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Iwasaki Wf] and DBW2 results exhibit a clean exponential fall-off for all values 
of Ng where data are available. This implies that the effect of exponentially localized 
near-zero modes vanishes within numerical accuracy. In these cases the density of 
near-zero modes should be very small, and their average localization length should 
be of order one in lattice units; hence one is well outside the Aoki phase. 

The "hamiltonian" H{—M, 05) also serves as a kernel of the generalized overlap 
operator obtained in the ^ 00 hmit, and A'^ and I are the two quantities which 
control the locality of this operator. The small but non-zero value of m^es, found using 
the Wilson gauge action at quenched (3 = 6.0 when both A^^ and the separation are 
large, suggests that in this case the localization scale of the limiting operator might not 
be sufficiently small compared to the physical scale of the simulation. In contrast, 
the Iwasaki and DBW2 results imply a highly localized effective four-dimensional 
operator. 

When overlap fermions are employed, as with domain-wall fermions it is neces- 
sary to verify whether locality of the overlap operator is obtained, and whether the 
localization scale is small compared to the physical length scale. This should become 
a routine practice in any overlap simulation! 

The picture for the phase diagram of QCD with Wilson fermions we painted 
in this paper is for a good part conjectural. In particular, we did not prove the 
hypothesis about the mobility edge presented in Sect. 4.2. In addition, some of our 
more rigorous arguments are based on the assumption that, in the regions where the 
mobility edge does not vanish, a maximum localization length l^nax exists. However, 
we would like to emphasize again that our conjecture appears to be the simplest 
possible way in which we can understand the collected evidence about the quenched 
and unquenched phase diagrams, incorporating both numerical [13 CHI 1201 HI] and 
analytical P 123 CIl IH 123 results. 

In order to test our proposal for the phase diagram, it would be interesting to 
study the Ward identity (|2.9p numerically, in particular in the quenched two-flavor 
theory. The second term on the left-hand side should show the 1/mi divergence 
already in finite volume. A comparison with the pionic condensate will test whether 
the l/rrti divergence saturates the Ward identity for very small momenta. A study 
of the dependence of all terms in the Ward identity on the volume should make it 
possible to see in which region of the phase diagram Goldstone excitations occur in the 
infinite- volume limit. Obviously, such a numerical study will have to include a twisted 
mass 1211 ; (exceptional) configurations should be discarded. Numerical studies 

can be extended to study localization lengths of near-zero modes as well as non-zero 
modes outside the Aoki phase, once the location of the various regions in Fig. 1 has 
been established. The localization length of non-zero modes would yield information 
on the value of the mobility edge. (As explained earlier, similar information may 
be obtained from (quenched) domain- wall fermion simulations.) Similar studies can 
also be done in the unquenched theory, but in that case the spectral density of the 
localized modes will depend in a more complicated manner on the randomness of 
the gauge field because of the back-reaction through the fermion determinant. It is 
interesting that the quenched theory provides a conceptually simpler arena to test 
the validity of our conjecture! 
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Appendix A 

In this Appendix, we consider the issue of localization from a somewhat differ- 
ent angle. Instead of using our criterion for localization of individual eigenmodes, 
Eq. 1)4.11) . we use Anderson's criterion of (partial) absence of diffusion jSZj. The 
argument reviewed here is originally due to Ref. [IH], see also Ref. j2H|- 

In the language of condensed matter, the hermitian Wilson-Dirac operator H = 
H{mo) is a "tight-binding" hamiltonian. This means that the "electrons" can reside 
only on the sites of a regular lattice. This hamiltonian lives in four space dimen- 
sions, with coordinates x or y, and determines the evolution of the system in a (fifth, 
continuous) time dimension, with coordinate s. We will assume that the electron 
encounters a random potential on each lattice site, defined by some probability dis- 
tribution. Since there is no back reaction of the electrons on the random potential, 
this situation corresponds to the quenched approximation of QCD, where of course 
the gauge fields provide the random potential. We introduce the advanced Green 
function Q{x, y; s) and its Fourier transform Q{x, y; A) 

g{x,y;s) = e{s)J2^^{x)^i{y)e-'^"\ 

n 

G{x,y-\) = ^v]/„(a;)vi/t(y)(A-A„)-S (A.l) 

n 

which solve the equations 

{id/ds-H)g{x,y;s) = i5\x - y)5{s) , 

{X-H)g{x,y;X) = 6\x - y) . (A.2) 

To carry out the Fourier transform, A has to be in the upper half plane, whereas 
taking it in the lower half plane corresponds to Q{x,y; A) being the transform of the 
retarded Green function. 

We now adopt the point of view that the electron is partially localized if there 
is a non-zero probability for the electron to be found at a specific location y after 
an arbitrarily long time, given that it started out at some other location x. This 
probability is given by {\Q{x,y; s)\'^) for s — oo, where (. . .) denotes the statistical 
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average over the random potential at each lattice site (or, in the QCD case, at each 
lattice link). The idea is that this non-zero probability arises from localized states, 
whereas the probability to escape to infinity comes from extended states. Obviously, 
in a finite volume there is no clear distinction between these two probabilities, but 
one may study the dependence of {\Q{x,y; for s — oo on the volume. If it stays 
non-zero in the infinite- volume limit, we have (partial) localization. 

Following Ref. we may express the limiting value of {\Q{x,y] s)p) in terms 
of ensemble averages of the Fourier transform Q{x,y] X): 

\im{\g{x,y;s)\') = \im n T ds e-^^ {\g{x,y; s)\') (A.3) 

s— >oo \ / rj^O JO ^ ' 

= \im/j2^^{x)^i{y)^,{y)^l{x)——^ — \ (A.4) 

= lim ^ [dX {g{x, y;X + zr//2) g{y, x; X - zr//2)) . (A.5) 

Equality ()A.3|) follows after a partial integration (after which the limit rj ^ can be 
taken); the other two follow from elementary integrations. These equations are valid 
in finite volume, and we will assume that they are valid in infinite volume as well. 

We observe that the integrand on the right-hand side of Eq. ()A.5|) corresponds 
precisely to Eq. ()3.1|1 . i.e. the pion two-point function, for A = and identify- 
ing 77/2 = mi. With the Anderson criterion for localization, we immediately con- 
clude that the A- integral of the two-point function has a. l/rj divergence provided 
lim^^oo {\g{x,y', -s)!^) > 0. In finite volume this is always the case, whereas in infinite 
volume, this is the case if the electron is (partially) localized. 

Anderson localization implies that, by definition, the probability for the electron 
not to escape to infinity, given by lim^^oo {\g(yX,y; s)p), is non-zero, and thus that 

{g{x,y;X + ir])g{y,x;X-ir]))(xr]'^, 1] (A. 6) 

for a range of A (note that this correlation function is positive). It follows that 
localization for a given A may be defined as the existence of this divergence [IHl- This 
coincides with our analysis based on the mobility-edge hypothesis: if the mobility 
edge Ac > 0, all modes with A < Ac contribute to the 1/mi (or l/r]) divergence, and 
thus the right-hand side of Eq. ()A.5|1 does not vanish, leading to a non-zero probability 
for the electron not to diffuse. 

In reality, it may not be easy to measure the limiting value of {\g{x,y; s)\'^) 
experimentally, because all the electronic states are filled up to the Fermi energy. 
The nature of the eigenstates close to the Fermi energy determine the macroscopic 
properties (at sufficiently low temperature). A disordered system is an insulator (zero 
electric conductivity) if the electronic states at the Fermi energy are localized, whereas 
a metallic behavior (non-zero conductivity) is observed if the electronic states at the 
Fermi energy are extended. A metal-insulator phase transition occurs when a mobility 
edge reaches the Fermi energy. The Aoki phase transition is in this sense a special 
kind of a metal-insulator transition. For a discussion employing these concepts in the 
context of the chiral phase transition in continuum QCD, see Ref. [HUj . 
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Finally, we mention that the concept of an effective lagrangian for the long-range 
degrees of freedom which is used to study the phase diagram of lattice QCD [T2[I221, 
is also widely used in condensed-matter disordered systems, see Refs. 123 ISO] and 
references therein. 



Appendix B 

Using the operator formalism and the notation of Ref. the numerator in 
Eq. ()6.11|) is a sum of four positive terms. For a fixed gauge field (and in an arbitrary 
normalization) one of these terms is 

(O' |q(0, 0) f cl(x, r) ct(x, r) f ^^/^ c|(0, 0)| O') . (B.l) 

Here T is the second- quantized transfer matrix, and c are fermion creation and 
annihilation operators. The up and down arrows represent two different flavors. Spin 
indices have been suppressed. The state |0') is a reference state annihilated by all 
the c's. It encodes the boundary conditions in the fifth dimension. (For notational 
simplicity we have given the result when the quark mass is zero; all arguments gen- 
eralize to the case of a non-zero mass, as well as to the other three terms that we do 
not show.) 

The composite operator c| (x, r) C| (x, r) belongs to the pseudo-scalar density 
J^qix, t) located in the middle of the five- dimensional bulk. Therefore, to be reached 
from one of the boundaries, one has to act Ns/2 times with the transfer matrix 
(we assume that Ng is even). If an admissibility condition is imposed, the (first- 
quantized) transfer matrix T{M,a^) has a gap [i.e. there are no eigenvalues in some 
open neighborhood of one). In the limit A^^ ^ oo^ j'Ns/2 i^g^omes proportional to the 
ground-state projector \0h) {Oh\- (The ground state \0h) of T is obtained by filling 
the Dirac sea of states that correspond to all eigenvalues a;„ > 1 of T(M, as).) In 
this case expression ()B.1|) is proportional to (0^1 c|(a;, r) C|(x, r) |0i^), which is in fact 
identically zero, because the (in general non-zero) states c-^ \0h) and cj^ \0h) have dif- 
ferent flavors and, hence, are orthogonal. Thus, the numerator in Eq. (jfi.lll) vanishes 
and chiral symmetry is recovered (exponentially) in the limit A^s — ^ oo. The case of 
realistic gauge actions is discussed in Sect. 6. 
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